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Abstract 

We present preliminary measurements of branching fractions for the charmless two-body decays 
B° — > 7r + 7r _ and K + tt~ , and a search for £?° — > K + K~ using a data sample of approximately 227 
million BB decays. Signal yields are extracted with a multi-dimensional maximum likelihood fit, 
and the efficiency is corrected for the effects of final-state radiation. We find the charge-averaged 
branching fractions (in units of 10 -6 ): 

B(B°^tt + tt-) = 5.5 ±0.4 ±0.3, (1) 
B(B° -> K+TT-) = 19.2 ±0.6 ±0.6, (2) 
B(B° K + K~) = <0.40. (3) 

The errors are statistical followed by systematic, and the upper limit on K + K~ represents a 
confidence level of 
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1 INTRODUCTION 



Charmless hadronic two-body B decays to pions and kaons provide a wealth of information on CP 
violation in the B system, including all three angles of the unitarity triangle. The time-dependent 
CP asymmetries in the tttt system can be used to measure the angle a ^Q; the decay rates for 
the Kir channels provide information on 7 [2J; and the time-dependent CP asymmetry in ir°Kg 
approximately measures (3 in the standard model |3j and is a sensitive probe of new physics in 
the b — > s penguin-decay process 4 . Recently, direct CP violation in decay was established in 
the B system through observation of a significant rate asymmetry between B — > K + ir~~ and 
B° —> K~tt + p>)@. As B-physics experiments accumulate much larger data sets, charmless two- 
body B decays will continue to play a critical role in testing the standard model description of CP 
violation. 

In order to extract the maximum information from these decays it is necessary to understand 
the underlying hadron dynamics, and measurements of branching fractions for all of the charmless 
two-body B decays involving combinations of it^ 1 , K^, 7r°, and K® are invaluable in testing the 
various theoretical approaches [7] . We present preliminary measurements of branching fractions for 
the decays (Hj B° — ► 7r + 7r~ and K + 7r~, and a search for the decay B° — > K + K~ using a data set 
2.5 times larger than the one used for our previous measurements of these quantities Tabled 
summarizes previous experimental measurements [SJ 1101 111) and current theoretical estimates of 
the branching fractions for these decays. 

2 THE BABAR DETECTOR AND DATASET 

The data sample used for this search contains (226.6 ± 2.5) x 10 6 T(45) — > BB decays collected by 
the BABAR detector ^2] at the SLAC PEP-II e + e~ asymmetric-energy storage ring. The primary 
detector components used in the analysis are a charged-particle tracking system consisting of a 
five-layer silicon vertex detector and a 40-layer drift chamber surrounded by a 1.5-T solenoidal 
magnet, and a dedicated particle-identification system consisting of a detector of internally reflected 
Cherenkov light (DIRC) providing K~n separation over the range of laboratory momentum relevant 
for this analysis (1.5-4.5 GeV/c). 

Table 1: Summary of existing branching fraction measurements (in units of 10 -6 ) and theoretical 
estimates for the decays B° — > 7r + 7r~, K + tt~, K + K~ . Theory estimates are from Beneke et al. 
and Keum in Ref. 0. 



Mode 


B (BABAR) 9 


fi(Belle) 10^ 


£(CLEO) 11 


Theory 




4.7 ±0.6 ±0.2 


4.4 ±0.6 ±0.3 


+ 1.4 +0.5 
*'°-1.2 -0.4 


4.6-11.0 


K + ir- 


17.9 ±0.9 ±0.7 


18.5 ±1.0 ±0.7 




12.7-21.0 


K+K- 


< 0.6 


< 0.7 


< 0.8 


0.007-0.080 
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Table 2: Summary of total detection efficiencies (%) for signal decays determined in GEANT Monte 
Carlo samples without FSR effects, compared with the results using PHOTOS and the leading- 
order QED calculation. We use the latter result in calculating the branching fraction and take the 
difference with PHOTOS as the systematic uncertainty. Uncertainties are statistical only. 



The data sample used in this analysis is similar to that used in the BABAR measurements of direct 
CP violation in K + ir~ 5 and time-dependent CP-violating asymmetry amplitudes S W7r and ^3] 
(the reader is referred to those references for further details on the analysis technique). Relative 
to the event selection applied in the CP analyses, we remove the requirement on the difference in 
the decay times (At) between the two B mesons in order to minimize systematic uncertainty on 
the branching fraction measurements. All other selection criteria are identical to those applied in 
Refs. (HIE]- We identify B — > h + h~ (h = tt or K) candidates with selection requirements on track 
and Cherenkov-angle (0 C ) quality, l?-decay kinematics, and event topology, and determine signal 
and background yields through a multi-dimensional maximum-likelihood fit. The final sample 
contains 69264 events and is defined by requirements on the energy difference, \AE\ < 150 MeV, 
and energy-substituted mass, 5.20 < tties < 5.29GeV/c 2 , of the selected B candidates [T4"] . 

The efficiency of the selection criteria is determined in large samples of GEANT-based Monte 
Carlo simulated signal decays. We include the effects of electromagnetic radiation from charged 
particles using the PHOTOS simulation package [15 \. The addition of final-state radiation (FSR) 
leads to the development of a low-energy tail in the distribution of AE for B° — > h + h~ signal 
candidates, which can cause some fraction of events to fail the \AE\ < 150 MeV requirement. We 
have implemented a detailed QED calculation |16| up to 0(a) in order to correct the efficiency 
obtained by the PHOTOS simulation. Table [2] summarizes the comparison of the efficiencies for 
the different modes assuming no FSR, the PHOTOS result, and the QED calculation. For the 
branching fraction measurement we use the efficiency as determined by the QED calculation, and 
take the difference with respect to PHOTOS as the systematic uncertainty. 

In addition to signal 7r + 7r~, K + tt~ , and (possibly) K + K~ events, the selected sample includes 
background from the process e + e~ — > qq (q = u, d, s, c). Possible backgrounds from other B decays 
are small relative to the signal yields (< 1%), and are treated as a systematic uncertainty. We use 
an unbinned, extended maximum-likelihood fit to extract simultaneously signal and background 
yields in the three topologies (tttt, Kit, and KK). The fit uses the discriminating variables rrtESi 
AE, 8 C , and the Fisher discriminant T described in Ref. [01, where the likelihood for event j is 
obtained by summing the product of the event yield and probability V% over the signal and 
background hypotheses i. The total likelihood for the sample is 



Mode No FSR PHOTOS QED 



7T+7T- 40.9 ±0.2 39.9 ±0.2 39.4 ±0.2 
K+tt- 39.9 ±0.2 38.9 ±0.2 38.4 ±0.2 
K + K- 38.6 ±0.3 37.8 ±0.3 37.6 ± 0.3 



3 ANALYSIS METHOD 
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The probabilities V{ are evaluated as the product of the probability density functions (PDFs) with 
parameters dTj, for each of the independent variables Xj = {tues, AE, T, 0~}, where 6^ and 
6~ are the Cherenkov angles for the positively- and negatively-charged tracks, respectively. The 
largest correlation between the Xj is 13% for the pair (wies, AE) and we have confirmed that it 
has negligible effect on the fitted yields. For both signal and background, the K^tt^ yields are 
parameterized as n K ± 7T ^ = iiKn (1 T Ak-k) /2, and we fit directly for the total yield uku and the 
asymmetry Ak-k- The result for Ak-k is used only as a consistency check and does not supersede 
our previously published result ■ 

The eight parameters describing the background shapes for rriES; AE, and T are all allowed to 
vary freely in the maximum-likelihood fit. We use a threshold function ^Fj\ for thes (1 parameter), 
a second-order polynomial for AE (2 parameters), and a sum of two Gaussian distributions for T 
(5 parameters). For the signal shape in m,Es, we use a single Gaussian distribution to describe all 
three channels and allow the mean and width to vary freely in the fit. For AE, we use the sum 
of two Gaussian distributions (core + tail) , where the core parameters are common to all channels 
and allowed to vary freely, and the tail parameters are determined separately for each channel from 
Monte Carlo simulation and fixed in the fit. Given that the tail is dominated by FSR effects, we 
take the shape directly from the Monte Carlo samples after correcting for the difference between 
PHOTOS and the QED calculation. For the signal shape in we use an asymmetric Gaussian 
function with different widths below and above the mean. All three parameters are determined in 
Monte Carlo simulation and fixed in the maximum-likelihood fit. The 9 C PDFs are obtained from a 
sample of approximately 430000 D* + — > D°tt + (D° — > K~~tt + ) decays reconstructed in data, where 
K T /tt^ tracks are identified through the charge correlation with the tt^ from the D* ± decay. The 
PDFs are constructed separately for K + , K~ , tt + , and ir~ tracks as a function of momentum and 
polar angle using the measured and expected values of 6 C , and its uncertainty. We use the same 
PDFs for signal and background events. 

Table 01 summarizes the fitted signal and background yields, and Kir charge asymmetries. We 
find a value of Ak-k consistent with our previously published result, and a background asymmetry 
consistent with zero. The signal yields are somewhat higher than the values reported in Ref. [2] 
due to the removal of the At selection requirement and the addition of the radiative tail in the 
signal AE PDF. In order to quantify the effect of FSR on the fitted yields, we perform a second fit 
using a single Gaussian for the AE PDF allowing the mean and width to vary freely. The results 
are shown in the second column of Table 03 where we find that ignoring FSR lowers the tttx yield 
by 4.5% and the Ktt yield by 2.4%. 

As a crosscheck, in Figs. ^ and 13 we compare the PDF shapes (solid curves) to the data using 
the event-weighting technique described in Ref. For each plot, we perform a fit excluding the 
variable being plotted and use the fitted yields and covariance matrix to determine a weight that 
each event is either signal or background. The distribution is normalized to the yield for the given 
component and can be compared directly to the assumed PDF shape. For ities, AE, and T , we 
find excellent agreement for signal nir and Ktt events (Fig. as well as the sum of all channels 
for background events (Fig. [5J). We have verified separately that the background PDF shapes agree 
for all three channels. Figure 01 shows the likelihood ratio Cs/J2^i f° r all 69264 events in the 
fitted sample, where Cs is the likelihood for a given signal hypothesis, and the summation in the 
denominator is over all signal and background components in the fit. We find satisfactory agreement 
between data (points with error bars) and the distributions obtained by directly generating events 
from the PDFs (histograms). 



10 



Table 3: Summary of the branching fraction fit using a sample of approximately 227 million BB 
pairs. For comparison, we show the results using a single Gaussian for the signal AE PDF, which 
would correspond to an analysis that ignores FSR effects. 



Parameter 


Nominal Fit 


Ignoring FSR 




491 ± 35 


469 ± 34 


Nk-k 


1674 ± 53 


1634 ± 52 




-0.135 ±0.030 


-0.135 ±0.030 


Nrk 


3.0 ±13.1 


5.3 ± 12.6 




32977 ± 194 


32998 ± 194 


NbKn 


20761 ± 169 


20801 ± 169 


AbK-K 


0.002 ± 0.008 


0.002 ± 0.008 


NbKK 


13358 ± 126 


13356 ± 126 



BaBar 

Preliminary 




5.27 5.275 5.28 5.285 5.29 



m ES (GeV/c ) 




5.285 5.29 



(GeV/c) 



| 200 



100 



1 1 

b) 

/ 


i 1 i 

[\BaBar - 

1- 1 Preliminary 


i , 1 i * 


-0.1 


0.1 
A E (GeV) 




0.1 
A E (GeV) 



50- 




-2 2 

Fisher Discriminant 




-2 2 

Fisher Discriminant 



Figure 1: Distributions (points with error bars) of mES) AE, and J~ for signal 7r + 7r - (a,b,c) and 
K + ir~ (d,e,f) decays using the weighting technique described in Ref. [TSJ. Solid curves represent 
the corresponding PDFs used in the fit. The distribution of AE for K + ir~ events is shifted due to 
the assignment of the pion mass for all tracks. 



11 




Figure 2: Distributions of a) tues, b) AE, and c) T for qq background events (points with error 
bars) using the weighting technique described in Ref. ^H]. Solid curves represent the corresponding 
PDFs used in the fit. 




0.5 1 0.5 1 0.5 1 

kk Likelihood Ratio K71 Likelihood Ratio KK Likelihood Ratio 



Figure 3: The likelihood ratio Cs/J2 A> where Cs is the likelihood for each event to be a signal irir 
(left), Ktt (middle), or KK (right) event. The points with error bars show the distribution obtained 
on the fitted data sample, while the histograms show the distributions obtained by generating signal 
(red) and background (blue) events directly from the PDFs. 
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Table 4: Summary of relative systematic uncertainties on yields, efficiencies, and number of BB 
pairs. For the K + K~ yield we show the absolute uncertainty. The total uncertainties for 7r + 7r~ 
and K + n~ are calculated as the sum in quadrature of the individual contributions. 



Source 


7T+7T 


K + tt- 


K+K- 


yields 


3.8% 


1.8% 


6.8 


efficiency 


2.6% 


2.5% 


2.0% 


V- 


1.1% 


1.1% 


1.1% 


Total 


4.7% 


3.3% 


n/a 



4 SYSTEMATIC STUDIES 

Systematic uncertainties on the branching fractions arise from uncertainties on the selection ef- 
ficiency, signal yield, and number of BB events in the sample. Uncertainty on the efficiency is 
dominated by track reconstruction (1.6%) and the effect of FSR (1.3%), which is taken to be the 
difference between the efficiency as determined in the PHOTOS simulation and the QED calcula- 
tion (Table |2j). Uncertainty on the fitted signal yields is dominated by the shape of the signal PDF 
for T (2.9% for vrvr, 1.5% for Kit) and potential bias (2.2% for vrvr, 0.9% for Kit) in the fitting 
technique determined from large samples of Monte Carlo signal events and a large ensemble of 
pseudo-experiments generated from the PDF shapes. Uncertainties due to imperfect knowledge of 
the PDF shapes for tties, AE, and C are all less than 1%. Table 21 summarizes the total uncer- 
tainty on the branching fractions, which is calculated as the sum in quadrature of the individual 
uncertainties. 

5 RESULTS and SUMMARY 

Table El summarizes the preliminary results for the charge- aver aged branching fractions. For com- 
parison, we use the efficiencies and signal yields determined under the assumption of no FSR and 
find B(B° -> 7T+71-) = 5.1 x 10" 6 and B(B° K + tt~) = 18.1 x 10~ 6 , which are consistent with 
our previously published results Taking into account FSR effects leads to an increase of the 
branching fractions by approximately 8% and 6% for tttt and Ktt, respectively. The upper limit on 
the signal yield for KK is given by the value of iVo for which £ max dN/ J °° £ max dN = 0.90, 
corresponding to a one-sided 90% confidence interval. Here, £ max is the likelihood as a function of 
N, maximized with respect to the remaining fit parameters. We find Nq = 25.9, and the branch- 
ing fraction is calculated by increasing the signal yield upper limit and reducing the efficiency 
by their respective total errors (Table 0J) . For the purpose of combining with measurements by 
other experiments, we have also evaluated the central value for the branching fraction and find 
B(B° K+K-) = (4 ± 15 ± 8) x lO" 8 . 

In summary, we have presented preliminary updated measurements of charge-averaged branch- 
ing fractions for the decays B° — > 7r + 7r~ and B° — > K + it~ , where FSR effects have been taken into 
account. We find a value of Axn consistent with the result in Ref. 0, and branching fractions 
6-8% higher due to the effect of FSR on the efficiency and signal-yield determination. This differ- 
ence should be taken into account when comparing with previous measurements of these quantities 
(Table ^) that do not include these effects. Our results are consistent with current theoretical 
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Table 5: Summary of branching fraction results in a sample of (226.6 ±1.2) x 10 6 BB pairs. We 
show signal yields Ng, total detection efficiencies (e) and branching fractions B in units of 10 -6 . The 
errors are statistical and systematic, respectively, and the upper limit on B° — > K + K~ corresponds 
to the 90% confidence level. 



Mode 


N s 


e(%) 


£(10- B ) 


7T + 7T~ 


491 ±35 ±11 


39.4 ±0.2 ±0.9 


5.5 ±0.4 ±0.3 


K + it- 


1674 ± 53 ± 15 


38.4 ±0.2 ±0.8 


19.2 ±0.6 ±0.6 


K+K- 


3.0 ± 13.1 ±6.8 (< 25.9) 


37.6 ±0.3 ±0.8 


< 0.40 (90% C.L.) 



estimates using various techniques [7]. We find no evidence for the decay B° — > K + K and set an 
upper limit of 4.0 x 10 -7 at the 90% confidence level. 
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